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SUMMARY 


The purpose of this report was the investigation of the symmetric modes 
of motion of an airplane and of the concepts which have been used to obtain 
relatively simpler expressions for the period of these motions. -The utility of 
the Reeves Electronic Analog Computer in this study was investigated. 

opecifically, the concepts concerning the symmetric motions were applied 
to the general body axes equations of motion in different ways. After the 
application of these simplifying assumptions, the modified equations were 
solved for the period of the motion concerned. Numerical data for a repre- 
sentative jet fighter type aircraft at 40,000 feet, Mach number .75, and at 
sea level, Mach number .25 was utilized in the solution by the different 
methods. These values for the period were compared to the true values as 
determined from the complete equations of motion. 

The possibility of increasing the accuracy of the period obtained in 
one of the methods was examined. 

A major objective of this report was the investigation of the utility 
of the analog computer in solving problems of this type. Mathematical treat- 
ment, while not difficult, is cumbersome and slow. The analog computer, on 
the other hand, after having been supplied with information, achieves a solu- 
tion almost instantaneously. It was desired to examine the effects of changed 


stability derivatives and coupling between the equations utilizing the computer. 





The conclusions reached as a result of this investigation are: 

1. The concepts concerning the short period symmetric motion of an air- 
plane when applied to the general equations provide an acceptable 
approximation to the period of that motion. 

2. Concepts concerning the period of the phugoid mode, when applied, 
provide accuracies ranging from 30 percent low to le percent high, 
depending upon the method of application. It was concluded, there- 
fore, that use of the existing simplifications does not accurately 
describe the period of the true motion because essential couplings 
are neglected. 

3. The Reeves Electronic Analog Computer provided rapid accurate solu- 
tions to many sets of differential equations which had to be solved. 
The utility of the computer in investigations such as this in unques- 
tionable. Variation of stability derivatives or coupling between 
equations could be accomplished by changing one or several resistances 
and the new solution was immediately obtainable. 

This study was conducted at the University of Michigan, Ann Arbor, Michi- 

gan, during the period February to June 1955 under the supervision of Professors 


M. A. Brull and J. D. Schetzer of the Aeronautical Engineering Department. 





INTRODUCTION 


The advent of flexible high speed aircraft has somewhat complicated the 
problems which must be solved by the aerodynamicist and the structures spec- 
ialist. With these complications has come an increasing tendency to write 
the equations of motion of the airplane in the body axes system of coordi- 
nates. 

The purpose of this report was the investigation of the symmetric modes 
of motion of a representative modern airplane in the body axes system and in 
particular, the investigation of the validity of the simplifying assumptions 
which have been used to obtain rapid approximate solutions for the period of 
these motions. 

The assumptions which have been made are based upon the observations of 
Lanchester who was the first to call the long period motion the phugoid. 
Lanchester stated that phugoid motion takes place at essentially constant 
angle of attack and further that the pitching motion is so slow that the 
aerodynamic moment may be considered zero, or nearly so. When these assump- 
tions are applied to the wind axes equations of motion and additional terms 
Which do not materially affect the period are ignored, an expression is 
obtained for the period of the motion of .138 V. The same results have been 
obtained from considerations of kinetic and potential energy interchange + 

The current literature does not contain examples of the application of 


the phugoid simplifying assumptions in the body axes equations of motion. 











The short period mode of motion is defined as a motion in pitch of very 


short duration which takes place as the angle of attack,which has been changed 


due to a disturbancey returns to the equilibrium value. The motion is so rapid 


that it is assumed that the forward speed does not change and consequently, 


the motion is independent of changes in perturbations in the forward speed. 
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SYMBOLS 


aircraft velocity ft/sec 

aerodynamic force in X body axis direction lbs 

aerodynamic force in Z body axis direction lbs 

aerodynamic moment about the Y body axis ft lbs- 

perturbation velocity in X body axis direction ft/sec 

perturbation velocity in Z body axis direction ft/sec 

perturbation velocity in Y body axis rad/sec 

dot over a symbol - derivative of symbol with respect to time 
@perator fer derivative with respeey to tine, = 

dimensional derivative of the aerodynamic force in the X body axis 
direction with respect to the perturbation velocity u; i sec/ft 


OX ap sec/ft 


Ow 


OZ 
Se lb sec/ft 


oz 


als sec/ft 
Ow 


92 Ib sec/ft 

= ft 1b sec/ft 
OM Ft 1b sec/ft 
OM Ft Ib sec/ft 


OM rt di sec/radian 





angle of attack of M.A.C. relative to wind - degrees 

angle of pitch of X axis measured from the horizon - degrees 
steady state linear velocity along X axis - ft/sec 

steady state linear velocity along Z axis - ft/sec 

mass in slugs 

moment of inertia Aete Y body axis - slug ft 

Mach number 


Lift coefficient 


Sum of forces in X body axis direction 





SYMMETRIC EQUATIONS OF MOTION 


(BODY - AXES SYSTEM) 


The six equations of motion break down into two sets of three equations. 
One set describes asymmetric motion. The other set describes symmetric motion. 
The linearized equations for symmetric motion, neglecting Xq> in body axes 
coordinates are: 
= [X,-md]u + [X,]w + | <i. COSTES q 


[Zy]u+ [Z.-m]w + [ZqtmJ,-BB S22 Si)q 


[MjJu + [MyMyNw + [My - ly » Ja 


= QO 


My 


I 


0 


The determinant of the coefficients is set equal to zero and the stabi- 
lity quartic is obtained 


AXN* + BAS + CAX2 + DA + B= O 





where 
Az=l 
2% Mg | Ms 
=e © Ge oee 
M W M.. 
aig (tha) + ary ONeutla) + Te € sin Ga 
— 
+ (MXM) + BE (Kaze Zale) 
My 
D = Seay (2akw-Kuee) - SE (yeanMalan) + GES (MaZa-Maee 





is aa 
: ks 7 =a. a =L (Myth, 


B= BE EL (Meda Mate) + SEED (hi -Ma) 


The equations of motion are developed in Ref. l. 





When the analog computer is to be utilized, it is convenient to rewrite 
the equations of motion as follows: 


Xyu + Xyw - m Wiq-Mg cos a1 | 


9b be — 
my = Zu + Zw + [Zg+mU |}, - mg sin a, | 
ea = My + Myw + Myw + Mgq 


The computer circuit for the solution of these equations is included as 
Fig. 1. The airplane numerical data supplied to the circuit is controlled 
by positioning the potentiometers. Each potentiometer controls a specific 


airplane parameter. These are; 


He a On Nees 1h Pike 
oO. xX ‘a mW 4 12, Z, 
Sia. al 8. mg sin 0, 15. 4g + mU 7 
ee 9. M, 14. Mg 
Ok Iyy 10. M, 15. Ve eos 83: 








APPLICATION OF CLASSICAL ASSUMPTIONS 
FOR SIMPLIFICATION OF THE SOLUTIONS 


TO THE MODES OF MOTION 


The roots of the stability quartic yield the characteristics of the 
short period and phugoid motions. Various simplifying assumptions con- 
cerning these motions have been introduced in the past. Accuracy obtain- 
able after the assumptions have been applied varies. Normally, the short 
period motion is so heavily damped and so rapid that it is of little con- 
sequence. It does become important if the period is within the normal 
human response lag and not heavily damped, primarily because a pilot may 
inadvertently reinforce the motion when he attempts to correct it. The 
short period mode is important when platform stability is a consideration, 
The phugoid motion is a long period slow oscillation with weak damping and 
of such little consequence that it may be negatively damped and not render 
an airplane unacceptable. Because of the apparent unimportance of these 
modes in piloted subsonic aircraft, approximations to their period and 
damping are justified and desirable. 

However, the importance of these modes of motion in pilotless missiles 
at very high speeds has not been completely investigated. When a relatively 
long slender missile encounters severe turbulence or a sharp gust, the pos- 
sibilities seem to be present for considerable control difficulties. These 
motions could intensify the difficulties. 

It follows then that it is important to have a method for the rapid, 


accurate determination of the symmetric modes. 








SHORT PERIOD METHOD 


The short period motion is obtained by assuming that the motion proceeds 
at constant forward velocity (u= 0). This follows if it is reasoned that by 
the time the forward velocity has changed, the motion has damped out and has 
ceased to exist. The assumption permits a reduction of the equations of 
motion to essentially a coupling between the perturbations w and q. The 


equations of motion become: 


F,= 0 = ([%,-m]w + [Zg+mU,- SB Sr)q 


M, = O = [MAM-A]w + [Mg - Iyy A]q 


Solution of these equations as written involves the solution of a cubic. 
The equation becomes a quadratic if the perturbation gravity force is assumed 
very small and ignored. The quadratic is: 
De een 0 
where: 


A= 


1 
no lt 
lyy m 





Zq+mU 
mlyy ( q i) 


C= SuMg Exe, (Z +mU; ) 
mlyy mlyy *° 4 


The utility of the analog computer is demonstrated by the fact that the 
computer circuit for the short period method is achieved by uncoupling poten- 
tiometers 11 and 12 in Fig. 1. This uncouples and removes from the circuit 
all contribution to the short period of the u perturbation velocity and the 
solution is obtained using the reduced Z force and moment equations. 


Ref. 2 shows a simplified short period method in the wind axes system. 








PHUGOID METHOD 


The long period phugoid motion is usually obtained by applying Lan- 
chester's original assumptions. Ref. 3 states that in a phugoid motion the 
angular pitching movement of the aircraft is at all times very slow and this 
implies that the aerodynamic pitching moment is nearly zero. Hence, the 
angle of attack is at all times very nearly that corresponding to zero 
pitching moment. This implies nearly constant angle of attack. Ref. 3 
further assumes that @ and C; is strictly constant. 

This investigation revealed that a choice is available for the method 
of enforcing zero change in angle of attack in the body axes equations. 

By definition: 


Q = arctan Lin 
7 


oa oa, 
doe Sw Ott a du, 


if dW, = perturbation velocity in Z direction = w 
and dU,-= perturbation velocity in X direction = u 
then 
NOD oo, + ey, 
OW, OU 
Noe Wcos @ -u sind 


Vp 


PHUGOID METHOD I 


The expression for the perturbation angle of attack is usually 


Simplified by assuming the equilibrium angle of attack is small, and 





because of this assumption, that 
Ao = 
Dp 
By this method, the DeRcarpe en @® is related directly to the pertur- 
bation velocity w. 
| When Lanchester's assumptions are applied to the equations of motion 
based on this near equality for AQ, and Ad = w = O, the stability quartic 
reduces to a quadratic which gives an approximate solution. The quadratic 
is obtained by ignoring the moment equation and deleting all terms in w in 


the force equations. Two equations remain in the variables u and q. These 


are: 
cos @ 
Fy = O= [X,-mJu  -  [nW) 8 =1)q 
= ot sin 9 
us =- O= [Z5] ut [Z,+mU,-“2 = sq 
The determinant of the coefficients set equal to zero yields: 
AXN= + BN + C = 0 
Where: 
A = [ZqtmU,] 
Z 
B= -[X, (= + U) + mg sin Q)+ 2Wa] 
C= g [X, sin 0,-4Z, cos 9] 


When the analog computer is utilized, the two equations are rewritten 
as follows; 


mi = X,u - mWiq - mg cos 04q 


(ZgtmU}4 = mg sin Giq - Zu 
The computer circuit is shown in Fig. 2@. Specific airplane parameters 


are controlled by the potentiometers as follows: 


aden mM 5 2. mg sin O41 
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re 8. mg cos 9, 
9. mW 
Switch in Fig. 2 is open. 
This method of simplification produces a period which is greater than 


the true period. 
PHUGOID METHOD II 


A second method for the reduction for the complete equations to obtain 
the period of the phugoid mode was developed by the writers during this in- 
vestigation. This method involves the application of the identical simpli- 


fying assumptions of method I, but the complete equality 


Aq = Wcos@-usin@ io 
ts 
was used to enforce no change in angle of attack. Then w = u tan @ was 


substituted in the equations and the third equation was ignored. The three 


equations reduce to two equations which are: 


FL = 0 = [X,-m\ + X, tan @ju - [m Wy +B OS Pia 
F =0O06c [Z,+Z,,tand-mAtand Ju zs [Zq+mU, BB SAE4]q 
The determinant of the coefficients set equal to zero yields: 
AN + BA + C = 0 
where: 
Z, 
A =m [U,+W, tan a + 3] 


- [ (Zeta )( <2 =H tan a) + mg (sin ©, - cos 6, tan @) 


td 
I 


+ Wy (Z,+Z, tan a) ] 


aba 





Cc = g [sin 9, (X, tana@+ X,) - cos 0, (AZ, tan a)] 
The equations are rewritten as follows when the analog computer is used: 
mi = (Xy+Xy tan @) u - m Wiq - mg cos 0, q 
(Zgtm U1)q = mg sin 6) q - (ZytZy tan Q) u + mi tan a 
The computer circuit is shown in Fig. 2. Specific airplane parameters 


are controlled by the potentiometers as follows: 


1. Zg + my Ce ee Cane 
2. mg sin @, (- Xy + Xy tan a 
3. m 8. mg cos @, 

4, m tan @ 9. mW, 


De as 

Switch in Fig. 2 is closed. 

This method of simplification produces a period which is less than the 
true period. 

Solution by Method II provides a period which is exactly equal to .138 
Vp, arrived at in wind axes in Ref. 2, and is exactly equal to 4.44 2 de- 
rived in Ref. 3. Thus Method II is the true application of Lanchester's 
Basic Assumptions in body axes. These assumptions are based on the pure 
interchange of kinetic and potential energy. 

Ref. 4 states, "To secure a lift proportional to the square of the speed, 
the wings must exert a constant lift coefficient, so that their incidence must 
be constant. This implies that the aeroplane is statically stable as regards 
pitching, that it has a negligible moment of inertia about the pitching axis, 


and the length between wings and tail is small compared to the radii of cur- 


vature of the flight path." 


a 





Method II does not substantiate this statement since substitution of 
w= u tan ad, (for the condition da = 0) in the moment equation causes My 
and M, to be grouped in the same coefficient and these terms exactly cancel 
each other because they become equal in magnitude and opposite in sign. 

( [Ma | = IM, tan eye This leaves the moment equation without a static sta- 
bility contribution and produces the same effect as zero static stability. 
In fact, solutions of pairs of the simplified equations with the moment 
equation included as one of the pair show the motion to be aperiodic. This 
is the motion to be expected from an airplane with zero static stability. 
The condition of zero static stability may be solved quickly on the analog 
computer by positioning the two potentiometers which control M, and M, on 
the total airplane circuit to zero. Computer solutions for this case show 
aperiodic motion. 

It is apparent that the Lanchester idealization enforces neglect of the 
moment equation because as far as the phugoid mode is concerned, the moment 
equation becomes inconsistent with the two force equations. This requires 
that the two force equations be used to obtain the solution for the phugoid 


mode and that an error in period must be accepted. 


MODIFIED PHUGOID METHOD I1 


It should be noted that Aq #& a = Q, if applied in the body axes 
p 


equations, produces a period which is larger than the true period; whereas, 

Method II produces a period which is smaller than the true period. 
Comparison of the determinantal coefficients resulting from both Ad 

approximations was made and the terms in each of the coefficients which 


largely effected the period were examined. It was noted that the gravity 
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contribution in the caqefficient multiplying \ to the first power (Method IT) 
was completely destroyed by tan @. 
[mg (sin @, - cos 0, tan a)] a= @, 

The classical explanation of phugoid motion is based on the interchange 
of kinetic and potential energy. The complete removal of all effects of 
gravity, therefore, seems illogical. In short, if the factor tan @ is re- 
moved from the term mA tan @ (second equation), the gravity effects are 
restored and the solution for the period of the phugoid is accurate to with- 
in one percent. The resultant period is nearly the actual value. 

This result is an experimental result arrived at on physical grounds 
alone and perhaps may not be applied with such accuracy in general because 
it cannot be mathematically justified. As shown above, the idealizations 
of Lanchester as generally applied are also not mathematically justified. 
(When A@ = O, the airplane is not simultaneously statically stable in the 
equations of motion.) If the additional reasoning concerning the gravity 
term is accepted as described on physical grounds, the coefficients of the 


Quadratic describing the phugoid mode are: 


A = m [Uz4+W, + 24 j 

B= - [(Zqtni ) (A tan a) + mg (sin 9, - cos 0,) 
+ Wy (Z2)+Zy tan a) ] 

C = g [sin ©, (X,+Xy tan @) - cos 01 (ZytZy tan @)] 


The computer solution for Modified Method II is obtained by charging the 
setting of potentiometer 4 in Fig. 2.(switch closed) to remove the factor tan a. 

Numerical examples for Method I, Method II, and Modified Method II are 
included in this report. Method I is 12 percent high, Method II is seven 
percent low and Modified Method II is within one percent, all referred to 


the actual period. 
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PHUGOID METHOD III 


One other method for an approximate solution for the phugoid mode has 
been used. This method makes use of the fact that phugoid motion is essen- 
tially a coupling between u mr The method of simplification used is to 
delete the w coupling terms in the force and moment equations thus uncoupling 
the Z force equation from the system. The X force equation and the moment 


equation thus reduced are solved simultaneously. The equations are: 


Fy, = O= [X,-m]u - [mW +B COS Pi )q 
Mp =O= [M, Ju+l My - Ty > lq 


The resulting cubic equation is: 


AN? + DAZ + CN 4+D = O 
where: 

Az 

B [es 


m Iyy 


pe Mu g cos 2. | 
lyy 


This method produces an error in the period of 25 to 30 percent. 

Solutions on the Reeves Electronic Analog Computer show this error. 
The computer solution is obtained by deleting the connections to potentio- 
meters 6, 9, and 10 in Fig. 1. This effectively uncouples the Z force equa- 
iiLon, The most powerful term in reducing the period 
was determined to be [My]W. This was shown by step by step replacement of 
the deleted terms. When [M,]w alone was reconnected, the period was restored 


to the correct value. 


15 





One other idealization was applied to the equations of motion in an 
attempt to find a simple accurate representation of the phugoid mode. Because 
the phugoid is described as a long slow oscillation about a datum line, it 
Seemed reasonable to assume that motions of this type would be characterized 
by very small accelerations in the direction normal to the flight path. This 
may be expressed as: 

w-Uiq = 0 
This equation was solved for q and substitution for q was made in the equations 
of motion. Pairs of the equation were solved with varying results. Solution 
of the two force equations produced aperiodic motion; solution of the Z force 
equation and the moment equation produced a short period mode frequency; and 
solution of the X force equation and the moment equation produced a cubic. 
This cubic was not solved but known accurate values of the phugoid roots were 
substituted to ascertain whether or not the correct phugoid mode was contained 
mi the cubic. It was found that the roots did not satisiy the equation even 
remotely and the method was discarded. 

It seemed logical to assume further that if the normal accelerations are 
small or zero, that the condition of zero angle of attack change would occur 
simultaneously. When this additional condition was imposed, the equations of 
motion reduced to three equations in u but all the solutions indicated aperio- 


dic motion. 
NUMERICAL RESULTS AND DISCUSSION 


The following examples are selected to facilitate comparison of the 
results obtainable by the several methods discussed above. The airplane 
upon which this is based is a representative jet fighter type aircraft. 


Two case are presented: first, the case of Mach number 0.75 at 40,00 feet; 
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second, the case of Mach number 0.25 at sea level. Normal mathematical 
solutions are presented and then are verified by Reeves Analog Computer 
solutions for the first case, The second case will be presented as analog 


computer solutions only. 


CASE I 


AIRPLANE DATA AND VELOCITIES 


S = 600 sq ft Xy = -5.19 lb sec/ft 
m = 1180 slugs X,, = 48.8 lb sec/ft 
Ta = 80,000 slug ft2 Z, = -70.8 lb sec/ft 
h = 40,000 ft altitude Z.y = ~807 lb sec/ft 
a@= 2.48 degrees = 0, M, = -637 lb sec 
Vp = 728 ft/sec M: = 10.08 lb sec? 
U, = 728 ft/sec M, = 27.5 1b sec 
Wy = 31.5 ft/sec Zq =-4.94 x 10% lb sec 
M = 0.75 Mg = -128.4 x 103 lb ft sec 


Using this data, the stability quartic coefficients are: 


As=l 

B = 2.2054 
C = 6.8699 
D = .O4604 
E = .02291 


By Bairstow's approximate factorization and one iteration for the phugoid 
roots, the quadratic is? 

N+ 00564 + 100332 = O 
The period of the phugoid is 109 seconds and the time to damp to half 


amplitude is 246 seconds. 
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The short period roots were obtained from the remaining factor of the 

quartic: 

N= + 2.19976. + 6.85419 = O 
The period of the short period mode is 2.65 seconds and time to damp to half 
amplitude is 0.63 seconds. 

The complete solution was obtained by using the Reeves Electronic Ana- 
log Computer and these results are shown in Fig. 4. The values obtained by 
the computer method and the stability quartic method are seen to be in excel- 
lent agreement and these values will be used as the basis for comparison of 


the various methods employed for simplification. 
SHORT PERIOD METHOD 


It is assumed that the short period mode is essentially a coupling be- 


tween the perturbations w and q alone. By this method, the equation is; 


AN? + BXN+C = O 
where: 
ee aa 
B= faneo 
Gowler 


From this, the period is 2.66 seconds and the time to damp to half amplitude 
is 0.63 seconds. This is in excellent agreement with the results obtained 


from the quartic. The computer solution is shown in Fig. }. 
PHUGOID METHOD I 


The second method to be compared is the previously discussed phugoid 


simplification in which AQ was considered to be equal to y° The quadratic 
Pp 
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equation for this method is: 


AN2 + BO+C = O 
where: 
A = 8.5306 
B= .04336 
C= (02275 


The period by this method is 122 seconds and the time to one-half amplitude 
is 273 seconds. The period is seen to be about 12 percent in error. The 


computer solution is shown in Fig. 5. 


PHUGOID METHOD II 


The third method to be compared is the method which requires that when 


da = 0, w=utanad, The quadratic equation which results is: 


AM +BX+C = O 
mere : 
A=. 6 156 
B= .0555 
C= aoo 


The method gives a period of 101 seconds and the time to damp to half ampli- 
tude is 213 seconds. The computer solution is shown in Fig. 6 for further 
comparison. It is seen that by this method the period is about eight percent 


too low. 


MODIFIED PHUGOID METHOD II 


The fourth method to be compared is the method evolved in this report. 


This method requires that when A@= 0, w=utanad. The quadratic equation 


7 





1 els 


AM= +BXY+C = O 
where: 
A= 8.92 
Hs 
C= 034 


The period by this method is 110 seconds and the time to damp to half ampli- 
tude is 29.6 seconds. The computer solution is shown in Fig. 7 for further 
comparison. It is seen that by this method the period is accurate to less 


than one percent. 
PHUGOID METHOD IIT 


The fifth and final numerical solution is the last phugoid case where 
w is considered to be zero and the second force equation is ignored. This 


method produces the cubic equation: 


AN? + BAe 4 Cha Dee eo 
where: 
A=1 
B = 1.60939 
C = .01788 
C = .01106 


This produces one negative real root and a complex pair. The period by this 
method is 75.7 seconds and the time to damp to half amplitude is 204 seconds. 
The computer solution is shown in Fig. 8 for further comparison. The error 


by this method is about 30 percent. 
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CASE II 


AIRPLANE DATA AND VELOCITIES 


S = 600 sq ft X, = -6.3 lb sec/ft 
m = 1180 slugs Xy = 165 lb sec/ft 
Iyy = 80,000 slug ft2 Za, = -162.6 lb sec/ft 
h = sea level Zy = -980 1b sec/ft 
Q@= 6.5 degrees = 0, M, = -1345 lb sec 
Vp = 279 ft/sec M: = -102.7 lb sec* 
Uz = 279 ft/sec M, = 153.7 Ib sec 
Wi = 31.9 ft/sec Z4 = -4325 lb sec 
M= .25 Mg = -115000 lb ft sec 


Computer solutions for the second case, Mach 0.25 at sea level are 


shown in Figs. 9 through 143and the results are seen to parallel Case I. 
CONCLUSIONS 


Comparison of the numerical results and computer solutions of the 

several methods presented lead to the following conclusions: 

1. The concepts concerning the short period symmetric motion of an 
airplane when applied to the general equations provide an acceptable 
approximation to the period of that motion. 

2. Concepts concerning the period of the phugoid mode, when applied, 
provide accuracies ranging from 30 percent low to l2 percent high, 
depending upon the method of application. It was concluded, there- 
fore, that use of the existing simplifications does not accurately 


describe the period of the true motion because essential couplings 


er 





are neglected. 

The Reeves Electronic Analog Computer provided rapid accurate solu- 
tions to many sets of differential equations which had to be solved. 
The utility of the computer in investigations such as this is unques- 
tionable. Variation of stability derivatives or coupling between 
equations could be accomplished by changing one or several resistances 


and the new solution was immediately obtainable. 
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FIG ~ 3 
CompuTER TRACES - COMPLETE EQUATIONS — /VaCH.75 
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